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1.  Introduction 


A  qutrit  density  matrix  is  of  order  3  and  depends  on  8  parameters.  There  are  many 
attempts  at  visualizing  the  state  space  of  a  qutrit1-4: 

1)  Using  all  8  Gell-Mann  matrices,  which  form  a  complete  set  for  expressing 
3x3  SU(3)  matrices; 

2)  Using  6  Gell-Mann  matrices  supplemented  by  different  matrices  in  place 
of  2  diagonal  ones;  and 

3)  Bloch  matrices  and  their  principal  minors. 

The  Bloch  matrix  approach  leads  to  many  complicated  constraints  on  the  qutrit 
parameters.  Previous  work  in  this  area  focused  on  understanding  the  structure  of 
8-D  parameter  spaces  through  2-  and  3-sections.  On  the  other  hand,  the  3-D  Bloch 
ball  state  space  of  qubits  are  easily  visualized.  Developing  geometrical  tools  to 
visualize  the  state  space  of  qutrits  in  a  similar  way  will  be  very  useful.  In  this  report, 
we  present  a  scheme  for  representing  qutrit  state  space  (QtSS)  in  3-D. 

2.  Unitary  Transformations 


2.1  General  Qutrit  States 

The  most  general  qutrit  density  matrix  can  be  represented  using  the  special  unitary 
group  SU(3)  invariant  form  as 

p  =  i/3+n.l  (1) 

Here  /3  is  3  x  3  unit  matrix,  n  =  (m,  m,. . ns)  are  8  real  parameters,  and  A  =  (At, 
2.2,. . .,  A  s)  are  3  x  3  Gell-Mann  matrices.  Using  them  one  obtains 


i  +  n3  +  l 

ni  ~  i-n2 

n4 

-  in5 

+  in2 

3  3  yf3 

n6 

—  in7 

(2) 

n4  -1-  ins 

i 

2  n8 

n6  +  in7 

3 

“  W  - 

2.2  Matrix  Exponential  for  Density  Matrix 

Given  a  Hamiltonian  H  and  an  angle-like  parameter  9,  the  corresponding  unitary 
transformation  is  given  as  U  —  el9p.  In  general  p  is  a  Hermitian  matrix.  For 
traceless  and  Hermitian  H  matrices  the  exponentials  are  widely  known.  The 


Approved  for  public  release;  distribution  is  unlimited. 

1 


calculation  becomes  more  involved  for  a  general  3x3  matrix,  which  is  the  case 
with  the  density  matrix  under  consideration. 

The  characteristic  equation  (CE)  of  p  derived  from  Eq.  2  is  given  as  follows: 

p3  -  p2  +  c2p  -  c3  =  0.  (3) 

Using  the  definition  r2  =  ^  ( a 2  +  qf),  the  coefficients  of  the  CE  are  given  by 

c2  =  ^{1  -  Tr(p2)}.  (4) 

c3  =  i{l  -  37r(p2)  +  27r(p3)}.  (5) 

Expressing  the  CE  in  terms  of  its  eigenvalues, 

ip  ~Po)(p  -P1XP-P2)  =  0-  (6) 

It  leads  to  alternate  expressions  for  the  CE  coefficients: 

c2  =  PoPx  +  P1P2  +  P2P0,  c3  =  PoPiP2 •  (7) 

So  the  eigenvalues  are  expressible  in  terms  of  the  density  matrix  elements.  These 
coefficients  are  further  constrained  by  the  following  inequalities: 

0  <  c2  <  i  0  <  c3  <  (8) 

One  can  derive  the  constraints  on  the  traces  also  by 

0  <  Tr(p2)  <  <  Tr(p3)  <  \  (9) 

It  is  possible  to  express  any  function  of  the  density  matrix  (including  the 
exponential  for  the  unitary  transformation)  using  the  method  of  projectors  as  given 
in  de  Zela.1  Let  pkand  \pk  >  denote  the  eigenvalues  and  eigenvectors  of  the  density 
matrix.  Denoting  the  3-D  identity  matrix  by/3,  we  can  write, 

h  =  1.1=0  I T-k  ><  Tkl  P  =  Yl=oPk\Pk  ><  PkV  (10) 

The  matrices  | pk  ><  pk\  ik  =  0,  1,  2  can  be  expressed  as  a  combination  of 
eigenvalues  (pk)  and  powers  of  density  matrix  ip).  Also  any  function  of  density 
matrix  can  then  be  written  as 

fiP)  =  IL  0  f(Pk)  I  Pk  ><  Pkl  (11) 
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(12) 


The  unitary  operator  takes  the  following  form: 

U  =  ewp  =Z2k=0eWpk\pkXpk\. 

2.3  Unitarily  Transformed  Density  Matrix 

The  density  matrix  transformed  by  the  unitary  operator  is  given  as 
UpU+  =  ewppe~wp  =  (£2k=0eWpk\ pk  ><  pk\)p(X2k=0e~wpk\pk  ><  /tfc|).(13) 

One  needs  to  use  a  relation  for  density  matrix  due  to  the  Cayley-Hamilton  Theorem. 

p3  =  p2  -  c2p  +  c3.  (14) 

Its  use  reduces  the  higher  powers  of  density  matrix  to,  at  most,  the  quadratic. 

2.4  Unitary  Transformation  for  Special  Cases 

Details  about  the  calculation  are  given  in  the  Appendix.  Here  we  enumerate  the 
results  under  a  common  scheme. 

Case  I:  Roots:  0,  0,  1 

CE:  P2(p~  1)  =  0 

Constraints  on  CE  coefficients:  c2  —  0,  c3  —  0 

Unitary  operator:  U  —  el9p  —  I3  +  (eld  —  l)p 

Transformed  density  matrix:  UpU+  —  el9p  pe~i9p  —  p 

Case  II:  Roots:  0,  p±  ,  p2  (p2  >  p±) 

CE:  p(p-  px){p  -p2)  =  0 

=  \  ~  a,  p2  =  l  +  «  (with  a  =  c2,) 

1 

Constraints  on  CE  coefficients:  c2  <  -  ,c3  =  0 

Unitary  operator:  U  —  elGp  —  el 9^2  |cos0cr  +  i  (~~~j  {p  ~ 

Transformed  density  matrix:  UpU+  —  ei0ppe^iep  —  p 

1 

Case  III:  Roots:  0,  px  =  P2=~ 

CE:  p(p  —  p-t)2  =  0 

1 

Constraints  on  CE  coefficients:  c2  =  c3  =  0 
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Unitary  operator: 

Transformed  density  matrix: 

Case  IV:  Roots:  p0  =  Mr  =  AC 
CE: 

Constraints  on  CE  coefficients: 
Examples: 

Unitary  operator: 

Transfonned  density  matrix: 

Case  V:  Roots:  p0,  p  1  —  M2 
CE 

Constraints  on  CE  coefficients: 
Examples: 

Also, 


U  —  eiQp  —  ei9/ 2 

UpU+  —  ei9p  pe~i9p  =  ew/2pe~w/2  =  p 

(M  -  Mo)3  =  0) 
cl  =  27c2 

Cc2,c3)  =  (U^7?).(U^).(Us7g).etc. 

y  _  ewp  _  ete^o 

UpU+  —  el9p0  pe~l9p0  —  p 

(M-Mo)(M-Mi  )2  =  0. 
c3=^[9c2-2  +  2(1-3c2)3/2] 


Mi  =  M2  =^(1 +  /?),/?  =  V1  _  3c2 


1— 4c2  +9Co 


^  =  3^"»  =  5(1-2«' 


Unitary  operator:  U  —  el0P  =  1  [el6,fil(p  —  p0)  +  el9p0(p1  —  p)] 

Transformed  density  matrix:  UpU+  =  el9ppe~l9p 


=  P  + 


(sinf$0/2'\ 

\JJ2~) 


c3 


g(2+/? 


/?2)P+-(1+/?)P2 


2.5  Unitary  Transformation  for  General  Case 

The  general  case  with  all  3  real  roots  (p0  <  px  <  p2)  leads  to  both  coefficients  of 
the  CE  being  nonzero.  So  we  get  a  complicated  expression  for  the  roots  of  the 
resulting  cubic  equation.  One  has  to  solve  the  following  equations  to  calculate  the 
density  matrix  and  related  functions. 

h  =  iMo  ><  Mol  +  iMi  ><  Mil  +  IM2  ><  M2L 
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(15) 


and 


P  =  PolPo  ><  Pol  +Pi\Pi  ><  Pi\  +  P2IP2  ><  P2I 
P2  =  PolPo  ><  Pol  +  M1IM1  ><  Pil  +  P2IP2  ><  P2I 


(16) 

(17) 


Their  solution  allows  one  to  write  the  unitary  transformation  operator  as 

U  =  ewp  =  ewp0\ p0  ><  ^0|  +  ei9pl\p1  ><  p±\  +  ewp2\ p2  ><  p2 \ 

_  eWn0  (Pi-P)(P2~P)  _  (po — P) (p2  —  P)  +  eWii2  (Po-P)(Pi-P) 

(Pl-Po)(P2-Po)  Ol-Po)(P2-Pl)  (P2-Po)(P2-Pl)’ 


Rewriting  Eq.  14, 


(Mi  -  M0XM2  -  Mo)(M2  ~  Pi)U  =  el9p0a  -  el9^b  +  el9p2c. 

(19) 

with 

a  —  (m2  Mi)(M2  -  M)(Mi  -  p), 

(20) 

b  =  (p2~  M0XM2  -  p)(Mo  “  P). 

(21) 

and 

c  =  (Mi  -  Mo)(Mi  -  p)(Mo  “  p)- 

(22) 

Then  the  transformed  density  matrix  can  be  expressed  as 

(Mi  -  Mo)2(M2  -  Mo)2(M2  -  Pi)2UpU+ 

—  -(-  ei9p2c][e~i9p0a  —  e~i9pib  +  e~i9p2c\p 

=  a2  +  b2  +  c2  —  2abcos (ji1  —  p0)8  —  2bccos(p2  —  p±)6  +  2cacos(p2  —  Po)6 
=  (a  -  b  +  c)2  +  Aabsin2  {pl~Po)d  +  4hcsm2  ^2~Pl)e  -  Acasin2  (M2~Mo)g.  (23) 

We  use  the  following  relations, 

Mo  +  Mi  +  M2  =  I.  c2  =  P0P1  +  M1M2  +  p2Po>  c3  =  P0P1P2 ,  (24) 


to  find  that  each  of  the  combinations  ah,  be,  and  ca  vanishes.  One  also  gets 
(Pi  —  PoYiPz  ~  Mo)2 (M2  —  Mi)2  =  (a  —  b  +  c)2.  Finally,  the  transformed  density 
matrix  reduces  to  the  following: 

UpU+  =  e^Ppe-iQp  =  (25) 


This  rather  surprising  result  says  that  the  density  matrix  p  with  3  real  and  distinct 
eigenvalues  describes  a  pure  state. 
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3.  Effect  of  Unitary  Transformations  on  3-D  Vectors  for  Qutrit 
States 

For  a  density  matrix  to  represent  a  physical  qutrit,  it  must  satisfy  some  constraints: 

1)  Positivity  of  diagonal  elements 

°<i±n3+|<l,0<i-^<l.  (26) 


2)  Normalization 


Tr{p )  —  1. 


3)  Constraint  on  the  length  of  the  state  vector 

n.n  <  -. 

3 


4)  Nonnegativity  of  the  determinant 

detp  >  0. 


(27) 


(28) 


(29) 


Earlier  attempts  at  visualizing  these  general  constraints  focused  on  8-D  Cartesian 
space  as  shown  by  Sarbicki  and  Bengtsson,2  Mendas,3  Goyal  et  al., 4  and  Bengtsson 
et  al.5  The  2-  and  3-sections  of  these  8-D  objects  (only  2  and  3  nonzero  elements) 
have  been  found,  classified,  and  characterized  by  these  authors.  Beyond  3 
dimensions  it  is  impossible  to  visualize  the  state  space,  so  attempts  have  been  made 
to  find  alternative  schemes.  In  this  work,  such  a  scheme  is  presented. 


3.1  Qutrit  Constraints  for  Spin-1  Representation  as  3-D  Vectors 

Following  the  formalism  of  Kurzynski  et  al.6  and  Kurzynski7,  an  alternative 
representation  of  the  qutrit  density  matrix  based  on  the  symmetric  part  of  the  2  qubit 
is  our  starting  point.  Define: 

di  =  (<7j  +  icq)/ 2,  dt  =  (tq  -  icq)/2  .  (30) 


Then  the  Bloch  matrix  representing  a  spin-1  state  is  given  by 


COi 

d3 

d2 

—d3 

0)2 

d i 

~d2 

d-^ 

0)3 

(31) 


The  parameters  in  this  representation  are  connected  to  spin-1  observables. 

=<  Sf  >=  Tr(pSf), 
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(32) 


at  =<  Si  >=  Tr(pSi),  (33) 

and  qk  =<  Sfi  +  SjSi  >=  Tr{p(SiSj  +  5)5))},  k  *  i,j.  (34) 

Comparison  of  this  parametrization  of  this  qutrit  density  matrix  with  the  earlier  one 
based  on  the  Gell-Mann  matrices  gives  the  following  identities: 

%  =  ^  <73.^2  =  \a3>nA  =  i q2,n5  =  i a2,n6  =  -\qltn7  =  ^alt  (35) 


and  n3  =i(<u1-ft)2),n8  =V3  -i  +  i(<u1  +  ft)2}].  (36) 

Then  the  constraints  in  Eq.  20  take  the  following  forms: 

1)  Positivity  of  diagonal  elements  implies 

>  0,  (jl)2  >  0,  (jl)3  >  0  .  (37) 

2)  Trip)  =  1  implies 

uq  +  (jl>2  +  <u3  =  1.  (38) 

3)  Tr  (p2)  <  1  implies 

Z?=i(<u?  +  2r2)  <  1.  (39) 

Here  r f  =  -  (a2  +  qf)  and  equality  sign  hold  for  the  pure  states. 

4 

4)  det  p  >  0  implies 

[XLi(3oj?  -  2ojf  +  6cotr2)  -  A]  <  1.  (40) 


Here  A  =  -  ( a2a3q1  +  a3alq2  +  a1a2q3  —  q^qs),  and  equality  sign  hold  for 
the  pure  states.  The  relation  to  CE  coefficients  is  given  as 

c2  =  jU  -  Tr(p2)}  =  (<u2<u3  +  <u3<ui  +  "i"2)  -  E?=  irf,  (41) 

and  c3  =  ^{1  -  3Tr(p2)  +  2Tr(p3)}  =  -ZLrWjr^  +  A.  (42) 

It  is  apparent  that  the  3  constraints  can  be  rewritten  in  terms  of  the  following  3-D 
vectors: 

1)  u  =  (iq,  u2,  u3}  such  that 

u2  -  i{37r(p2)  -  1}  =  f  +  2r2)  -  \ 
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(43) 


iq  = 

J^K+2  +)-i 

(44) 

u2  — 

+2r22)-i 

(45) 

and 

u3  — 

+  2  r,2)-i 

(46) 

2)  v  =  {rq,  v2,  v3}  such  that 
v2  =  ^{9Tr(p3)  —  1}  = 

|Zi=1{<yf  +  3ri2  -  3coiri  +  3A}  -  i 

(47) 

>1 

+  3(1  -  ao-^r2  +  A}  - 

(48) 

v>  =  4 

+  3(1  -  co2)r2  +  A}  -  — , 

(49) 

and  v3  = 

+  3(1  -  m3)r32  +  A}  -  — . 

(50) 

3)  w  =  {w1;  w,  w3]  such  that 

w2  —  co1  +  co2  +  co  3  =  1,  (51) 

and  w  -  {V^h,  V^3-  (52) 

These  vectors  have  the  desired  properties  such  that  for  pure  states  it2  =  1  =  v2, 
and  for  maximum  mixed  states  u2  =  0  =  v2,  and  they  also  obey  the  constraints 
0  <  it2  <  1  and  0  <  v2  <  1.  The  3rd  vector  w  is  always  of  unit  length.  The 
original  8  degrees  of  freedom  of  a  general  qutrit  state  is  distributed  as  3+3+2  among 
the  3  vectors.  They  are  always  in  the  positive  octant  of  a  3-D  sphere  of  unit  radius. 
The  mixed  states  reside  inside  the  volume  interior  to  the  “pure  qutrit  state  surface”. 
This  is  similar  to  the  Bloch  ball  of  a  qubit. 

3.2  Effect  of  Unitary  Transformation  on  the  QtSS  Vectors 

It  was  found  earlier  that  almost  all  the  cases  lead  to  no  change  in  the  density  matrix. 
The  special  case  with  one  double  root  with  CE:  (p  —  p0)(p  —  Mi)2  =  0  is  the  only 
case  when  this  does  not  hold  true.  The  transformed  density  matrix  is  given  as  the 
following: 

UpU+  =  ewPpe-wP  =  p  +  (fHwA  [c3-i(2  +  /?-/?2)p  +  i(l  +  /?)p2].  (53) 
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It  is  clear  from  the  definitions  that  B  —  u  and  c->  —  — —  -u2  +  —  v2,  and  so  the 

^  6  27  3  27 

connection  with  3-D  vectors  can  be  easily  established. 

4.  Conclusion 


In  this  work  we  have  shown  that  the  8-D  qutrit  state  space  can  be  alternatively 
visualized  in  3  dimensions  using  3-D  vectors,  and  unitary  transformations  affect 
their  lengths  and  directions.  In  the  future  it  will  be  interesting  to  see  if  the  QtSS 
3-D  vectors  offer  advantages  over  traditional  methods  for  visualizing  the  properties 
and  quantifying  measures  of  the  entangled  qutrits. 
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Appendix.  Unitary  Transformation  for  Special  Cases 
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Results  for  some  special  cases  are  as  follows: 

(I)  Both  coefficients  are  zero. 

c2  =  0,  C3  =  0.  (A-l) 

These  are  all  pure  states  with  eigenvalues  p  =  0,0,1. 

(II)  Second  coefficient  is  zero. 

c2  =£0,  c3  =  0.  (A-2) 

Which  gives  Tr(p3)  =  1  —  3 c2.  The  characteristic  equation  (CE)  for  this  case  is 
given  by 

P3  ~  M2  +  p2g  =  MO2  -  p  +  c2)  =  0.  (A-3) 

It  has  3  roots  out  of  which  one  is  always  zero.  The  other  2  roots  can  be  either  2  reals 
or  1  complex  conjugate  pair.  This  leads  to  the  following  special  cases: 

(Ila)  Two  real  nonzero  roots  (c2  <  1/4  or  Tr(p2)  >  1/2) 

Define  =  —  4 c2  ,  then  the  roots  are  0,  (1  +  a)/ 2,  and  (1  —  a)/ 2, 

and  p  =  diag(0,  (1  +  a)/2,  (1  —  a)/2).  (A-4) 

The  Cayley-Hamilton  Theorem  can  be  used  according  to  which  the  density  matrix 
itself  satisfies  the  CE.  Then  we  find  that 

U  =  eWp  —  ete/2[Cos(0a/2)  —  (i/a)sin(Qa/ 2)  +  (2/a)sin(6a/2)p].  (A-5) 
(lib)  One  double  root  (c2  =  1/4  or  Tr(p 2)  =  1/2) 

The  roots  are  a  =  0  and  p  —  diag{ 0, 1, 1).  There  is  a  single  root  of  0  and  a  double 
root  of  Vi.  The  unitary  transformation  matrix  is 

U  =  ewp  =7-2(1  -ei9'2)p.  (A-6) 

(He)  Pair  of  complex  conjugate  roots  (c2  >  1/4  or  Tr(p2)  <  1/2) 

Define  a'  =  ^/4c2  —  1  then  the  roots  are,  0,  (1  +  ia')/2,  and  (1  —  ia')/2, 
and,  p  =  diag( 0,  (1  +  ia')/2,  (1  —  ia')/ 2).  (A-7) 

We  also  get 

U  =  e10p  =  e10/2[cosh(0a'/2)  +  (i/a')sinh(0a'/2)  —  (2/ 
a')sinh(0a'/2)p].  (A-8) 
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(Ill)  Unitary  Transformation  for  case  with  2  nonzero  roots  (c2  <  1/4) 


The  general  unitary  transformation  is  given  by 


UpU+  =  ewppe-wP. 

(A-9) 

Here  p'  is  the  original  3x3  matrix  and  el9p  is  the  diagonal  matrix  exponential  for 
the  case  with  real  roots  (the  complex  root  situation  is  similar).  It  can  be  recalled 
that  for  a  qubit  this  results  in  the  rotation  of  the  3-dimensional  Bloch  vector. 
Rewriting 

aU  —  (Ax  +  SiP)  +  i(A2  +  B2p). 

(A- 10) 

Where 

A±  =  acos(9a/2)cos(0/2 )  +  sin(0a/2)sin(0/2), 

(A- 11) 

A2  =  acos(0a/2)sin(0/2)  —  sin(0a/2)cos(0/2), 

(A- 12) 

and 

=  —2sin(da/2)sin(8/2),B2  =  —  2sin(8a/2)cos(6/2). 

(A- 13) 

We  get 

a2UpU+  =  (A2  +  A\)p  +  (A1B1  +  A2B2){pp  +  pp)  + 
(Bl  +  B2)ppp  +  i(A1S2  -  A2B1)(pp  -  pp). 

(A- 14) 

with 

A\  +  A\  =  a 2  cos2(0a/2)  +  sin2(0a/2), 

(A- 15) 

B2  +  Bj  =  4sin2(0a/2), 

(A- 16) 

A1B1  +  A2B2  —  —  4sin2(0a/2), 

(A- 17) 

and 

A1B2  A2B1  —  2cc  sin  (y)  cos  (y). 

(A- 18) 

Also  let  D  —  (cos(0a/2)  +  (i/a)sin(9a/2)}.  Then  the  unitarily  transformed 
density  matrix  is  given  as 

D  e-Wa/z(-d3)  Deiea/2(-d^)' 

UpU+  =  DeWa'2(-d 3)  m2  eWa/2( -dj  .  (A- 19) 

De-iea'2{-d2 )  e-Wa/2(-dJ  co3 

Comparing  this  with  Eq.  31,  it  is  seen  that  under  a  unitary  transformation  the 
density  matrix  elements  get  multiplied  by  various  trigonometric  functions.  This 
induces  changes  in  the  lengths  and  angles  of  the  vectors  ( u ,  v,  w). 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


3-D 

3-dimensional 

8-D 

8-dimensional 

CE 

characteristic  equation 

QtSS 

qutrit  state  space 

SU 

special  unitary  group 
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